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Abstract 

The  effects  of  surface  polar  phonons  on  electronic  transport  properties  of  monolayer  graphene 
are  studied  by  using  a  Monte  Carlo  simulation.  Specifically,  the  low-held  electron  mobility  and 
saturation  velocity  are  examined  for  different  substrates  (SiC,  Si02,  and  HfCh)  in  comparison 
to  the  intrinsic  case.  While  the  results  show  that  the  low-held  mobility  can  be  substantially 
reduced  by  the  introduction  of  surface  polar  phonon  scattering,  corresponding  degradation  of  the 
saturation  velocity  is  not  observed  for  all  three  substrates  at  room  temperature.  It  is  also  found 
that  surface  polar  phonons  can  influence  graphene  electrical  resistivity  even  at  low  temperature, 
leading  potentially  to  inaccurate  estimation  of  the  acoustic  phonon  deformation  potential  constant. 

PACS  numbers:  72.80.Vp,72.10.Di,72.20.Ht,73.50.Dn 
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Since  the  experimental  realization  of  graphene  in  2004,-  its  extraordinary  electric  proper¬ 
ties  have  led  to  near  exponentially  growing  attnetion  in  possible  applications.  In  particular, 
the  linear  band  structure  in  the  vicinity  of  the  Dirac  points,  and  resultant  massless  fermions, 
have  excited  interest  in  both  its  applicability  for  electronic  devices  and  as  a  testbed  for 
quantum  held  theory  phenomena.—  The  extremely  high  intrinsic  mobilities  that  have  been 
observed  in  suspended  graphene  at  room  temperature,-  in  conjunction  with  its  true  two  di¬ 
mensional  structure  of  graphene,  naturally  lead  to  its  employment  in  higher  speed  graphene 
transistors  devoid  of  short  channel  effects.-  For  such  high  frequency  applications,  the  sat¬ 
uration  current  and  low-held  mobility  are  material  properties  that  need  to  be  accurately 
determined. 

When  a  graphene  layer  is  in  close  proximity  to  a  polar  substrate,  inelastic  carrier  scat¬ 
tering  with  surface  polar  phonons  (SPPs)  can  result  in  significant  reduction  in  the  low-held 
mobility  of  graphene.---  Due  to  the  inelastic  nature  of  SPP,  they  also  provide  pathway  to  cur¬ 
rent  saturation,  in  conjunction  with,  or  as  a  substitute  for,  intrinsic  optical  phonons.  While 
a  number  of  studies  generally  suggested  the  negative  influence  of  reduced  saturation  veloc¬ 
ities  (due  to  the  relative  small  energies  of  SPPs),--—  conflicting  reports  exist  that  predict 
a  very  different  picture  (including  enhanced  velocities)  based  on  the  analysis  of  Boltzmann 
transport  equation.-^ 

In  this  work,  we  investigate  the  influence  of  carrier  scattering  due  to  SPPs  on  the  elec¬ 
tronic  transport  properties  of  monolayer  graphene,  via  a  full-band  ensemble  Monte  Carlo 
method.  Particularly,  the  low-held  electron  mobility  and  saturation  velocity  are  calculated 
in  the  presence  of  three  different  substrates  (SiC,  SiC>2,  and  HfC^)  and  compared  with 
those  of  intrinsic  graphene  at  room  temperature.  Furthermore,  we  examine  the  impact  of 
SPPs  on  the  low  temperature  electrical  resistivity,  with  attention  to  its  implication  on  the 
experimental  determination  of  the  acoustic  phonon  deformation  potential  constant. 

The  Monte  Carlo  model  adopted  in  the  calculation  utilizes  the  complete  electron  and 
phonon  spectra  in  the  hrst  Brillouin  zone.  While  a  tight-binding  band  is  used  for  the  elec¬ 
tronic  energy  structure,  all  six  branches  of  the  graphene  phonon  spectra  are  considered  with 
the  phonon  dispersion  and  electron-phonon  scattering  rates  obtained  from  first  principles 
calculations.—  The  effect  of  degeneracy  is  accounted  for  by  the  rejection  technique,  after  final 
state  selection.—  The  distribution  function  is  obtained  self- consistently  from  the  ensemble 
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simulation.  The  SPP  scattering  rate  is  introduced  by  following^ 


27 r  *  ^ 


2jr2  e~2qd 


TS(ki)  h  ^  2 e(q)2  q 


(1  +  cos 9)(nq  +  -  ±  -)£(£/  -  Ei±  Tiujs)  , 


(1) 


where  q  =  |kf  —  kj|  is  the  SPP  momentum,  Ef  ( )  the  final  (initial)  electron  energies,  d 
the  distance  between  monolayer  graphene  and  the  substrate  (0.4  nm),  ujs  the  SPP  energy, 


and  E2  =  E^l 


re^+l  • 


the  square  of  Frohlich  coupling  constant.  The  (l  +  cos#)  factor 


originates  from  the  overlap  integral  of  the  pseudospin  part  of  the  electron  wavefunction  and 
e(q)  is  the  dielectric  function  in  graphene.  In  addition,  T  contains  the  dependence  on  high 
(low)  frequency  dielectric  constant  of  the  polar  substrate  Ed?  (k°s)  along  with  the  normalized 
area  A.  The  scattering  by  ionized  impurities  in  the  substrate  is  not  included  in  the  effort  to 
clearly  identify  the  role  of  SPPs.  The  values  for  relevant  substrate  parameters  can  be  found 
in  Refs.  6  s|. 

Figure  1  shows  the  SPP  scattering  rates  calculated  at  300  K  for  three  substrates,  SiC, 
Si02,  and  HfC>2.  For  comparison,  the  strength  of  electron-optical  phonon  interaction  inher¬ 
ent  in  graphene  (via  the  deformation  potential)  is  also  plotted  from  a  first-principles  analy¬ 


sis 
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The  results  clearly  illustrate  the  dominance  of  SPPs  over  intrinsic  optical  phonons  in 


graphene  for  the  entire  electron  energy  under  consideration.  Furthermore,  this  enhancement 
in  scattering  is  more  pronounced  for  the  substrates  with  small  ojs  that  is  apparent  from  the 
onset  energy  of  emission  process  (e.g.,  Hf02  -  19.4  meV;  Si02  -  60.0  meV;  SiC  -  116  meV). 
Due  to  the  Coulombic  nature,  the  SPP  scattering  is  a  function  of  electron  density  n  and  sub¬ 
sequent  screening  in  graphene.  Throughout  the  calculation,  we  assume  n  —  lx  1012  cm-2 
along  with  the  static  screening  function  e(q)  in  the  random-phase  approximation.—  The 
intrinsic  scattering  via  deformation  potential  interaction  has  no  dependence  on  n. 

Figure  2  shows  the  electron  drift  velocity  vs.  the  electric  field  at  T=300  K,  for  intrinsic 
graphene  and  graphene  on  different  substrates.  As  expected,  the  addition  of  electron-SPP 
interactions  lead  to  general  decrease  of  electron  drift  velocities  in  the  low-field  region.  Conse¬ 
quently,  the  low-field  mobilities  are  reduced  for  all  three  substrate  with  the  largest  decrease 
in  HfC>2.  However,  the  difference  with  the  intrinsic  case  becomes  progressively  smaller  with 
the  increasing  field  (thus,  the  increasing  average  electron  energy)  and,  in  the  case  of  SiC>2 
and  SiC,  the  velocity  appears  to  saturate  at  a  higher  value.  No  degradation  is  observed  even 
for  HfC>2.  Clearly,  this  behavior  does  not  follow  the  saturation  velocity  model  suggested 


by  Ref. 


9  based  on  the  onset  of  optical  phonon  emission,  intrinsic  phonon  or  otherwise 
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(i.e.,  SPPs).  Rather,  it  is  in  general  agreement  with  more  recent  studies  that  solved  the 
Boltzmann  transport  equation  numerically  assuming  a  displaced  Fermi-Dirac  distribution 
function 

The  apparent  inconsistency  of  increased  total  scattering  rate  and  absence  of  velocity 
degradation  in  the  high-held  region  may  be  explained  by  considering  the  linear  energy  dis¬ 
persion  and  the  characteristics  of  the  Frohlich  interactions.  As  it  is  well  known,  a  change 
in  electron  energy  (say,  via  an  inelastic  scattering)  does  not  automatically  relax  the  drift 
velocity  in  monolayer  graphene  near  the  K  and  K'  points.  What  matters  is  the  direction 
of  the  final  momentum.  Due  to  the  Coulombic  nature,  the  Frohlich  interactions,  including 
those  by  SPPs,  prefer  small  angle  events.  When  a  SPP  scattering  occurs,  the  electron  will 
likely  emit  (absorb)  a  SPP,  with  the  phonon  momentum  pointing  along  the  radial  direction 
toward  (away  from)  the  Dirac  point.  Consequently,  the  velocity  of  the  scattered  electron 
may  not  change  significantly.  In  contrast,  optical  phonon  scattering  of  intrinsic  graphene  is 
via  deformation  potential  interactions  and,  thus,  randomize/relax  effectively  the  direction 
of  the  electron  momentum/drift  velocity.  Ignoring  the  difference  between  the  Frohlich  and 
deformation  potential  interactions^  can  lead  to  inaccurate  depiction  of  transport  properties. 

When  the  electron-SPP  interactions  are  taken  into  account  along  with  other  scattering 
mechanisms,  they  provide  an  additional  channel  for  energy  and  momentum  relaxation.  Ac¬ 
cordingly,  the  average  electron  energy  becomes  substantially  lower  at  a  given  electric  field. 
As  the  applied  field  increases,  it  means  that  the  shift  in  the  distribution  function  (along 
the  direction  of  the  electric  field)  would  be  smaller  in  energy  along  with  a  shorter  tail  com¬ 
pared  to  that  in  intrinsic  graphene.  This  is  clearly  visible  in  Fig.  3  plotted  for  the  ky  —  0 
cross-section  at  20  kV/cm.  The  observed  shorter  tails  in  the  negative  kx  space  (i.e.,  with 
negative  velocities)  as  well  as  in  the  high  energy  region  with  nonlinear  dispersion^  appear 
to  more  than  compensate  the  additional  momentum  relaxation  of  SPP  scattering  on  the 
SiC>2  and  SiC  substrates,  leading  to  larger  saturation  velocities  shown  in  Fig.  2.  In  the  case 
of  HfC>2,  the  distribution  is  much  less  heated  with  a  significant  population  in  the  negative 
half  space.  Due  to  the  very  strong  inelastic  scattering  by  surface  phonons,  the  drift  velocity 
may  not  have  reached  the  saturation  point  at  20  kV/cm  and  could  grow  further.  Cooling 
of  electron  distributions  via  surface  phonon  interactions  may  also  reduce  the  self-heating  as 
the  generated  heat  (i.e.,  SPPs)  is  on  the  substrate  and,  thus,  can  be  more  readily  dissipated. 
Another  interesting  point  to  note  from  Fig.  3  is  that  the  electron  distribution  function  in 
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graphene  does  not  resemble  that  of  a  highly  degenerate  case,  certainly  not  when  subject 
to  an  appreciable  electric  held.  A  simple  approximation  of  a  displaced  box-like  functioned1 
cannot  describe  transport  properties  accurately. 

The  strong  influence  of  the  substrate  on  low-held  mobility  observed  at  300  K  raises  a 
possibility  that  the  SPP  scattering  may  be  an  efficient  mechanism  even  at  low  temperatures. 
This  is  a  distinct  possibility  due  to  the  small  phonon  energies  (Ug.  Figure  4  plots  the  electrical 
resistivity  p  in  graphene  vs.  temperature  for  all  four  cases  under  consideration.  Indeed,  the 
resistivity  is  substantially  affected  by  the  substrate  conditions  at  T  >  50  K  and  the  effect 
of  SPP  scattering  can  not  be  immediately  separated  from  those  of  other  mechanisms  (most 
importantly,  the  acoustic  phonon  scattering).  The  acoustic  deformation  potential  is  often 
determined  experimentally  utilizing  the  assumption  that  the  low  temperature  resistivity 
is  linearly  dependent  on  temperature  with  a  slope  that  is  proportional  to  the  square  of 
acoustic  deformation  potential.-—  However,  the  calculation  results  indicate  that  this  may 
not  be  possible;  the  extracted  deformation  potential  (Dac)  is  not  an  intrinsic  quantity  but  an 
’’effective”  parameter  that  includes  substrate  effects  and  screening  (i.e. ,  the  charge  density 
in  graphene)  among  others. 

Concerning  the  specific  results,  the  intrinsic  case  gives  Dac  &  6.8  eV  with  a  well-defined 
linear  region.  This  (i.e.,  the  strong  linearity)  is  due  to  the  negligible  contribution  of  optical 
phonon  scattering  at  low  temperatures  with  relatively  large  energies  (Tiuop  ~  160  meV). 
The  small  deviation  of  Dac  from  the  deformation  potential  constant  obtained  directly  from 
the  acoustic  phonon  scattering  rate  (~  4.5  eV)^  may  be  attributed  to  a  degenerate  electron 
density  considered  in  the  present  calculation  (n—  lx  1012  cm-2).  It  is  also  interesting  to  note 
that  this  value  (6.8  eV)  is  actually  rather  close  to  the  experimentally  extracted  on  a  non-polar 
substrate  (7.8  eV).1'  For  the  SiC  substrate,  we  estimate  Dac  «  7.1  eV  similar  to  the  intrinsic 
result.  As  T  >  150  K,  however,  there  is  a  substantial  increase  in  SPP  scattering  and  the 
difference  between  the  SiC  and  intrinsic  cases  becomes  more  discernible.  When  Si02  is  used 
as  the  substrate,  the  slope  of  p  is  further  increased  and  becomes  nonlinear  earlier  due  to  the 
small  SPP  energy.  The  deduced  value  in  the  linear  region  (T  125  K)  gives  the  appearance 
of  Dac  !=s  13.2  eV,  which  is  not  unlike  16—18  eV  estimated  experimentally  on  SiC>2.-  Finally, 
HfC>2  has  the  largest  effect  among  the  three  substrates.  The  resistivity  is  greatly  increased 
with  no  apparent  linear  region  in  the  temperature  range  under  consideration.  Accordingly, 
it  is  not  plausible  to  determine  Dac.  The  relevance  of  SPP  scattering  even  in  the  T  <  150  K 
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regime  can  explain,  at  least  in  part,  the  very  disparate  results  for  the  magnitude  of  acoustic 
phonon  deformation  potential  in  graphene. 5,8,13,17 

In  summary,  we  investigate  the  effects  of  substrate-induced  SPP  scattering  on  the  elec¬ 
tronic  transport  properties  in  monolayer  graphene  by  using  a  full-band  ensemble  Monte 
Carlo  simulation.  It  is  found  that  the  electron  velocity-field  characteristics  are  highly  de¬ 
pendent  on  the  choice  of  substrate  at  room  temperature.  Specifically,  the  Frohlich  nature 
of  the  interaction  appears  to  be  crucial  for  correctly  describing  the  saturation  of  drift  ve¬ 
locity.  The  simulation  also  shows  that  the  SPP  scattering  remains  an  efficient  mechanism 
even  at  low  temperatures  (T  >  50  K),  substantially  affecting  the  electrical  resistivity.  This 
makes  it  difficult  to  experimentally  determine  the  acoustic  deformation  potential  for  intrinsic 
graphene  in  close  proximity  to  a  substrate. 

This  work  was  supported,  in  part,  by  the  Defense  Advanced  Research  Projects 
Agency/HRL  Laboratories  (the  CERA  program),  the  US  Army  Research  Offce,  and  the 
SRC  Focus  Center  on  Functional  Engineered  Nano  Architectonics  (FENA).  MBN  wishes 
to  acknowledge  partial  support  from  the  Offce  of  Basic  Energy  Sciences,  US  DOE  at  Oak 
Ridge  National  Lab  under  contract  DE-AC05-00OR22725  with  LIT-Battclle,  LLC.  JMZ  ac¬ 
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Figure  Captions 


Figure  1.  (Color  online)  Surface  polar  phonon  scattering  rate  of  graphene  electron  on  SiC, 
Si02,  and  HfC>2,  for  the  electron  density  n  of  1  x  1012  cm-2  at  300  K.  Also  plotted  is 
the  intrinsic  graphene  optical  phonon  scattering  rate  at  300  K  (”  intrinsic” )  obtained  from 


Ref 
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Figure  2.  (Color  online)  Electron  drift  velocity  in  graphene  on  SiC,  SiC>2,  and  HfC>2.  The 
intrinsic  case  without  substrate  is  also  shown.  The  electron  density  is  1  x  1012  cnr2  at 
300  K. 


Figure  3.  (Color  online)  Electron  distribution  function  with  different  substrate  conditions: 
intrinsic  (dash-dotted),  SiC  (solid),  SiC>2  (dashed),  and  HfC>2  (dotted),  for  the  ky  —  0  cross- 
section  at  20  kV/crn.  The  box-like  function  corresponds  to  the  Fermi-Dirac  distribution 
displaced  by  the  Si02  SPP  energy  (60  meV)  in  a  simple,  metal- like  approximation  (i.e.,  EF 
kFT)  with  n  —  lx  1012  cnr2.  For  comparison,  the  equilibrium  Fermi-Dirac  distribution  at 
300  K  is  also  plotted. 

Figure  4.  (Color  online)  Electrical  resistivity  vs.  temperature  with  different  substrate 
conditions:  intrinsic  (circle),  SiC  (square),  Si02  (triangle),  and  Hf02  (diamond),  with 
n  —  lx  1012  cnr2.  The  slope  of  the  straight  lines  are  used  to  extract  the  acoustic  de¬ 
formation  potential. 
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FIG.  l:  Li  et  al. 
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FIG.  2:  Li  et  al. 
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FIG.  3:  Li  et  al. 
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FIG.  4:  Li  et  al. 
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